Mon. Not. R. Astron. Soc. 000, ITHT91 (2009) Printed 2 December 2009 (MN IMfeK style file v2.2) 



Cross-correlating WMAP5 with 1.5 million LRGs: a new 
test for the ISW effect 

U. Sawangwit 1 *, T. Shanks 1 !, R.D. Cannon 2 , S.M. Croom 3 , Nicholas P. Ross 1 ' 4 
: and D.A. Wake 1 ' 5 

1 Physics Department, University of Durham, South Road, Durham, DH1 3LE, UK 
2 Anglo -Australian Observatory, PO Box 296, Epping, NSW 1710, Australia 
3 School of Physics, University of Sydney, NSW 2006, Australia 

^Department of Astronomy and Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA 
' 5 Department of Astronomy, Yale University, CT 06520, USA 



Accepted 2009 November 16. Received 2009 November 06; in original form 2009 May 25 



ABSTRACT 

We present the cross-correlation of the density map of Luminous Red Galaxies (LRGs) 
and the temperature fluctuation in the Cosmic Microwave Background (CMB) as mea- 
sured by the five-year Wilkinson Microwave Anisotropy Probe ( WMAP) observations. 
The LRG samples were extracted from imaging data of the Sloan Digital Sky Survey 
(SDSS) Data Release 5 (DR5) based on two previous spectroscopic redshift surveys, 
the SDSS-LRG and the 2dF-SDSS LRG and QSO (2SLAQ) surveys designed to have 
average redshifts of z sa 0.35 and z ps 0.55. In addition we have added a higher-redshift 
photometric LRG sample based on the selection of the AAOmega LRG redshift sur- 
vey at z w 0.7. The total LRG sample thus comprises 1.5 million galaxies, sampling a 
redshift range of 0.2 < z < 0.9 over ps 7600 square degrees of the sky, probing a total 
cosmic volume of « 5.5/i -3 Gpc 3 . 

First, we find that the new LRG sample at z ps 0.7 shows very little positive 
evidence for the Integrated Sachs- Wolfe (ISW) effect. Indeed, the cross-correlation is 
negative out to ~ 1°. The standard ACDM model is rejected at ps2-3% significance by 
the new LRG data. We then analyse the previous samples at z ps 0.35 and z pa 0.55. 
As found by other authors, these results appear consistent with the standard ISW 
model, although the statistical significance remains marginal. W e also reprod u ce the 
same result for the magnitude limited SDSS galaxy samples of ICabre et al.l (|2006l) 
Taking the z ps 0.35 and z pa 0.55 LRG results in combination with the new z pa 0.7 
sample, the overall result is now more consistent with a null detection than with the 
standard ACDM model prediction. 

We then performed a new test on the robustness of the LRG ISW detections at 
z ps 0.35 and z pa 0.55. We made 8 rotations through 360° of the CMB maps with 
respect to the LRG samples around the galactic pole. We find that in both cases there 
are stronger effects at angles other than zero. This implies that the z pa 0.35 and 
z pa 0.55 ISW detections may still be subject to systematic errors which combined 
with the known sizeable statistical errors may leave the z pa 0.35 and z ~ 0.55 ISW 
detections looking unreliable. We have further made the rotation test on several other 
samples where ISW detections have been claimed and find that they also show peaks 
when rotated. We conclude that in the samples we have tested the ISW effect may 
be absent and we argue that this result may not be in contradiction with previous 
results. 

Key words: cosmology: observations - cosmic microwave background - large-scale 
structure of Universe. 



INTRODUCTION 



* E-mail: utane.sawangwit@durham.ac.uk 
f E-mail: torn. shanksOdurham. ac.uk 



Many observations now suggest that we l ive in a spatially 
flat, dark energy dominated Universe (e.g. iPerlmutter et al.l 



© 2009 RAS 



2 U. Sawangwit et al. 



1999 



2007 



Cole et al l 120051; iTegmark et all 120061 ; iRiess et all 



Dunklev et al.l 



2009) 



In such a cosmology, posi- 
tive correlation between the CMB and large-scale-structure 
(LSS) is expected due to the decaying gravitational potential 
jSachs fc Wolfe! 1 19671 ). The deviation of the CMB temper- 
ature in the vicinity of LSS is caused by the non-vanishing 
difference in the energy gained and lost by the CMB pho- 
tons as they traverse a region of over- or under-density. 
By integrating across all the potential wells along the line 
of sight from the surface of last scattering, the primodial 
fluctuations in the CMB are modified by this effect. This 
secondary anisotropy of the CMB is called the Integrated 
Sachs-Wolfe (ISW) effect and sometimes known as the late- 
time ISW effect since the dominance of dark energy in the 
cosmic energy budget at the present epoch is believed to be 
responsible for the current accelerating expansion and hence 
the decaying gravitational potential. For a spatially flat Uni- 
verse, a detection of the Integrated Sachs- Wolfe (ISW) effect 
would provide direct dynamical evidence of the accelerating 
expansion unlike the geometrical measurement inferred from 
standard candles such as the SNIa. 

The SNIa results coupled with the CMB evidence that 
the Universe is nearly flat, suggests there exists an exotic 
form of energy with negative pressure. The exact nature of 
this so-called dark energy is not yet known but it already 
entails many serious problems. Foremost amongst them are 
the fi ne-tuning problem and the cosmic co incidence problem 
(e.g. ICarrolill200ll ; I Peebles fc Ratrall2003l ). 

The ISW signal in the CMB-galaxy cross-correlation 
is very small, generally less than one jjK, and very 
difficult to detect. Previous ISW detections generally 
have less than 3a statistica l signifi c ance. These include 
the st udies of Fosalba et al. |2003l ). iPadmanabhan et alj 
i|2005l ) and ICabre et all (|200rj ) who used SDSS galaxies 
in both photo-z and magnitude limited samples and the 
WMAP3 dataset. Oth e r auth ors have used X-ray sources 
Boughn fc Crittenden! 120041 ') and NVSS radio 



Nolta et all 120041 '). Of these, it seems that up to now 



the most significant detection of the ISW effect comes 
from the NVSS radio so urces at 2.3cr. Other authors (eg 
iGiannantonio et al.l 120081 and iHo et all |200ST ) have made 
compilations of the other results and claimed up to 4a ISW 
detections in terms of the overall significance. The only other 
claims of ISW detections at high significance are the meth- 
ods that reduced the galaxy samples to focus only on regions 
of hig h or low underdensity. In particular, iGranett et al.l 
( 2008) cross-correlated the posit ions of » 100 superc lusters 
and voi ds in the MegaZ-LRG l|Collister et all 12007) sam- 
ple and iMcEwen et al.l l|2007l ) employed a similar wavelets 
method using radio sources from NVSS. 

Here we shall search for the ISW effect by using sam- 
ples of Luminous Red Galaxies (LRGs) from the SDSS DR5 
dataset. LRGs are the most luminous stellar systems in the 
Universe, residing in the most massive dark matter haloes. 
Having formed most of their stars much earlier and over 
short period of time, the objects appear red with reasonably 
uniform spectral energy distributions therefore these galaxy 
samples can be selected homogeneously and observed out to 
greater distance (or lookback time). Moreover, being mas- 
sive means tha t the LRGs are also a highly bias ed tracer 
of the LSS (e.g. lRoss et"al]|2007l . lWake et alj|2008l ). The se- 
lection techniques for z < 0.6 LRG samples have been well 



established in the literature. Many LSS studies have been 
carried out using these LRG samp les including the cl aimed 
detections of the ISW effect (e.g. ICabre et al.ll2006D . The 
recent spectroscopic survey bv lRoss et alj ( 20081 ) has shown 
that it is possible to extend the selection technique and 
hence the LRG sample out to z « 1. Applying this tested 
algorithm to the entire SDSS imaging significantly increases 
the effective volume and makes these LRGs ideal probes of 
large-scale structure. 

Our main goal is to detect the ISW signal in the CMB 
by cross-correlating WMAP5 map with the new z ~ 0.7 
LRG sample and to test the detection of the ISW effect 
caused by the LRGs at lower redshif t (z ~ 0.35, 0.55) as 
claimed by a number o f authors (e.g. IPadmanabhan et al.l 
120051 . ICabre et al.ll2006l ). These studies used the LRG can- 
didates extracted from the SDSS DR3 or DR4 whilst we are 
using DR5, ~ 50 percent and 20 percent increase in the 
area coverage, respectively. The larger sky coverage should 
provide a statistical advantage over the previous studies. 
Our new higher redshift LRG sample should also provide 
a chance to constrain the evolution if such an effect is in- 
deed detected. Moreover, a recent study by iDouspis et alj 
(2008) suggests that the ISW signal-to- noise can be opti- 
mised if the large-scale tracer probes out to a median red- 
shift of 0.8 but there is no further improvement after a red- 
shift of unity. The claim ap pears to be support e d by the 
cross-correlation analysis of IGiannantonio et alj (|2008l ) in 
which the signal-to-noise of the I SW detection from the 2 
Micron All Sky Survey (2MASS; I Jarrett et all |2000| ) is » 
4-6 tim es smaller than from the NRAO VLA Sky Survey 
(NVSS; ICondon et al.lll998T ) where z w 0.1 and 0.8 respec- 
tively, despite the fact that the two surveys have similar 
sky coverage and sky density (A^nvss ~ 2A^2mass). If this is 
true then our higher redshift LRG should be more sensitive 
to the ISW signal and will provide even higher significance 
of detection than previous studies using the LRGs which 
currently reach ~ 2a significance at best. The new sample 
therefore presents a fresh opportunity to test one of the most 
crucial manifestations of the accelerating expansion, obtain- 
ing independent confirmation of the geometrical inference 
of the SNIa result if detected and a challenge to the current 
standard picture of the Universe otherwise. 

The layout of this paper is as follows. We present the 
data in SjU We then outline the theoretical prediction and 
cross-correlation technique employed in this study in §3] and 
§3] respectively. The results and a range of analyses per- 
formed to ensure their robustness are given in §5] and §6] 
The additional sky rotation tests performed on our dataset 
and selections of previously claimed ISW detections are re- 
ported in SJ7] We then present the discussion and conclusion 
of our studies in Sj8] and §9] Throughout this study (unless 
otherwise stated), we assume a standard ACDM cosmology 
with J7a = 0.7, Sl m = 0.3, /baryon = 0.167, as = 0.8 and 
-ffo = lOO/i km s" 1 Mpc -1 (h — 0.7 where necessary). 



2 DATA 

2.1 CMB Temperature M.ap-WMAP5 

The CMB temperature maps used here are taken from 
the WMAP five-year data release (|Hinshaw et al.l 120091 ). 
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The data products are publicly available^] in Hierar- 
chical Equal Are a isoLatitude Pixelisation (HEALPix, 
iGorski et aT1l2005h format. Although the WMAP observes 
in five frequency bands, we choose to use only the three 
highest-frequency bands, namely, W at 94 GHz, V at 61 
GHz and Q at 41 GHz as the CMB anisotropy in these 
ranges are less susceptible to a contamination from the fore- 
ground anisotropy (i.e. synchrotron and free-free emission) 
than the lower frequency counterparts. This enable us to test 
for any wavelength dependence in the CMB-galaxy cross- 
correlation where one expects the ISW signal to be achro- 
matic. However, we shall concentrate our analysis mainly 
on the W band due to its relatively high resolution com- 
pared to the other bands, 12. '6 FWHM compared to 19/8 
for V and 29/4 for Q band. Despite the fact that the V 
band has lower noise than the W band (hence often the 
band of choice for this type of analysis), we do not observe 
any major difference in either the cross-correlation results 
or their statistical errors (see F ig. |3l). We also us e the Inter- 
nal Linear Combination (ILC. [Gold et all [20 09) to further 
check our results, although it should be noted that, accord- 
ing to the WMAP team, there could be a significant struc- 
ture in the bias cor rection map at scales smaller than « 10° 
|Limon et alj|2008h . 

We shall use the temperature maps at a resolution 
of Af s id e =512 (res=9) which for the whole sky, contains 
3 145 728 pixels each with an area of w 49 arcmin 2 . The 
foreground-contaminated regions of the sky, mainly in 
Galactic Plane and Magellanic Cloud including extragalac- 
tic point sources, are excluded using a co mbination of 'Ex - 
tended temperature analysis mask' (KQ75, Gold et ah! 20091 ) 
and 'Point source catalogue mask' IjWright et alj|2009 ). Af- 



Table 1. Summary of the LRG samples used in the cross- 
correlation analyses. 



ter applying the masks, we are left with 2 239 993 pixels 
(« 70 per cent). The maps contain thermodynamic tempera- 
tures with the dipole contribution subtract ed from the data 
by the WMAP team (|Hinshaw et al.ll2009h . 



2.2 Luminous Red Galaxies 

The Luminous Red Galaxy (LRG) photometric samples are 
extra cted from the SDSS DR5 (|Adelman-McCarthv et all 
12007ft imaging data based on three LRG spectroscopic red- 
shift surveys who s e med ia n redshifts are ~ 0.35 , 0.55 and 0.7 
jEisenstein et al.1 l200ll ; ICannon et"ai1 120061 ; iRoss et al.1 
l200ct ). In essence, these surveys utilised a crude but effec- 
tive determination of photometric redshift (photo-z), owing 
to the strong 4000 A break of a typical E/S0 galaxy spec- 
tral energy distribution (SED). As the break is redshifted 
through the SDSS u,g,r,i, and z bandpasses, its colour- 
colour track exhibits a distinctive turning point at various 
redshifts for different colour pairs. Moreover, their uniform 
SEDs ensure that they have an extremely tight locus in the 
colour space. This allows the potential LRGs in the desired 
redshift ranges to be selected uniformly using their loca- 
tions on the colour-colour plane coupled with the luminosity 
threshold set by the appropriate magnitude limit. 

These simple methods have been proven to be highly 
effective in selecting the intrinsically luminous early-type 
galaxies in the targeted redshift ranges as demonstrated 
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Figure 1. Redshift distribution of the three LRG samples in- 
ferred from the redshift surveys used in their selections. 



by t he SDSS-LRG, 2SLAQ and AAT-AAO redshift sur- 



veys l|Eisenstein et al.ll200ll ; ICannon et al.ll2006j ; I Ross et all 
2008). Although the LRG photo-z in these redshift ranges 
can be estimated quite accurately (|Padmanabhan et al.l 
120071 ; ICollister et al.ll2007h . we decided to base our study on 
the colour-magnitude cuts because a well defined photo-z 
error distribution is needed for the deconvolution to recover 
the real redshift distribution and could bias the analyses of 
the results. The colour-magnitude cut techniques used in 
the above spectroscopic surveys, applied to the entire SDSS 
DR5 dataset (only Northern Galactic Cap), results in ~ 1.5 
million LRG candidates and the redshift distribution of the 
survey is assumed for the corresponding photometric sam- 
ple. The outlines of the selection algorithms with the em- 
phasis on any differences in our criteria to that of the spec- 
troscopic surveys are given below (readers are referred to 
lEisenstein et al.ll200ll ; ICannon et aLll2006j ; IRoss et al.ll2008l 
for the detailed descriptions of the selection criteria). The 
number-redshift relations, N(z), (shown in Fig. [1]) used in 
the model predictions have been calibrated to include those 
differences. The summary of the three LRG samples is given 
in Table □ 

For the following sections, all magnitudes and colours 
are given in SDSS AB system (unless otherwise stated) and 
are cor rected for galactic extinction using t he galactic dust 
map of ISchleeel. Finkbeiner fc Davis! (|l998T ). 
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2.2.1 SDSS LRG 

The low redshift (median 2 w 0.35) LRG candidates 
are selected on the basis of their colours and magni- 
tude s following the ongoin g SDSS-LRG spectroscopic sur- 
vey l|Eisenstein et al]|200ll . E01 hereafter) which will con- 
tain more than 100 000 spectra and cover over lh~ 3 Gpc 3 
when completed. The survey is designed to be approximately 
volume-limited up to z « 0.4. The targets are selected us- 
ing g — r and r — i colour cuts with the magnitude limit, 
r pc tro < 19.5. Two sets of selection criteria (Cut I and Cut 
II in E01) are used to extract LRGs in two different (but 
slightly overlapped) regions of the gri colour space, sepa- 
rated by the turn-over of the gri colour track at z ~ 0.4. 

In addition to the criteria of E01, we also apply restric- 
tion on the bright limit in the r-band, i.e. r po t ro 17.5. This 
is mainly because Cut I is too permissive and allows under- 
luminous objects to enter the sample below redshift 0.2 and 
by imposing the bright limit, we restrict the sample to only 
2 > 0.2. The choice of r pc t ro 17.5 merits a brief explana- 
tion. The redshift-dependent luminosity threshold is imple- 
mented by one of the selection rules, r pc tro < 13.1 + cn/0.3 
(Eq. 4 in E01), where cy ~ g — r at z ~ 0.2 corresponds 
to g — r ~ 1.3 on the gri colour-colour track. This has 
been empirically confirmed to work sufficiently well using 
the spectroscopic sample, with only a few objects having 
r P ctro < 17.5 at z > 0.2 and vice versa. 

The LRG sample described above is then extracted from 
the SDSS DR5 imaging database using the SQL query by 
setting the flag PRIMTARGET to GALAXY_RED. This 
yields a catalogue of approximately 200 000 objects which 
after applying the bright flux cut mentioned above, be- 
comes 106 699 objects with the sky surface density of ~ 
13 deg" 2 . The average redshift of the LRG candidates as 
inferred from the spectroscopic sample of ~ 60 000 SDSS- 
LRG is z = 0.35 ± 0.06. 



2.2.2 2SLAQ LRG 

The 2dF-SDSS LRG and Quasar Survey (2SLAQ) is the 
spectroscopic follow-up of the intermediate redshift (z > 
0.4) LRGs from photometric data of the SDSS survey using 
the 2-degree Field (2dF) instrument on the 3.9-m Anglo- 
Australian Telescope (AAT). The survey was completed in 
2006 and contains approximately 13 000 spectroscopically 
confirmed LRGs with 0.4 < z < 0.8 in two equatorial 
strips covering w 180 deg 2 jCannon et alj feooe, and refer- 
ences therein). The primary and secondary sample of the 
survey (Sample 8 and 9, respectively) were selected using 
the SDSS g — r versus r ~ i colours in conjunction with the 
'de Vaucouleurs' i-band magnitude, (17.5 < id c v < 19.8). 
The colour cuts are similar to the Cut II used by E01 which 
targets the objects lie above the turning point of the early- 
type galaxy track in gri colour space. The turning point 
is caused by the 4000 A break moving into the r-band at 
2 w 0.4, making r — i colour increase rapidly whereas g — r 
remains nearly constant at 1.6-1.7 mag until z « 0.7. 

In order to qualify as 2SLAQ LRG candidates, objects 
are required to have d± ^ 0.65 for the primary sample 
and 0.55 ^ d± < 0.65 for the secondary sample, where 
d± = (r - i) - (g - r)/8.0. We shall only use the LRG 
candidates extracted following the primary sample cut, de- 



signed to target higher redshift candidates than the sec- 
ondary sample. We also utilise the star-galaxy separation 
criterion used by the 2SLAQ survey which has been proven 
to be very effective and the stellar contamination in the 
LRG sample is only 5 percent. The primary sample con- 
tains 67 percent of all 2SLAQ LRGs and has an average 
redshift of z — 0.55 ± 0.06. Applying the primary target 
selection, including the star-galaxy separation criteria, on 
the DR5 'best' imaging database in the NGC, a sample of 
655 775 photometrically classified LRGs is returned. Objects 
with BRIGHT or SATURATED or BLENDED but not DE- 
BLENDED flags are not included in our sample. 



2.2.3 AAQ, LRG 

Our new high redshift L RG sample is bas ed on the AAT- 
AAf2 LRG Pilot run (R oss et al. 2008, and references 
therein), using the 2dF instrument on the AAT. The sur- 
vey was carried out as a 'Proof of Concept' for a large LRG 
redshift survey, VST-AAf2 ATLAS survey. It was designed 
to target potential LRGs out to z m 1.0 with the average 
redshift of 0.7. Three different sets of selection criteria were 
employed in selecting the targets in order to test the AAS7 
spectrograph's ability to obtain reliable redshift with the 
minimum exposure time in average conditions. They ob- 
served over ss 10 deg 2 in three 2dF fields and the survey con- 
tains 1270 unique galaxy spectra with 804 high- confidence 
LRG redshifts. 

The selection rules used here follow the colour- 
magnitude cuts which utilise the riz colour plane and the 
'de Vaucouleurs' i-band magnitude. This selection forms the 
main part of the survey. In summary, the colour cuts exploit 
the upturn of the early-type galaxy track similar to that 
used by 2SLAQ and SDSS LRGs survey in selecting z > 0.4 
LRGs with gri colours. But in the riz colour space, the 
upturn occurs between redshift 0.6 and 0.7 as the 4000 A 
feature moves into the i band hence makes it ideal for select- 
ing potential LRG targets for the intended redshift range. 
The star-galaxy separation procedure uses the 2-band pho- 
tometry, akin to the method which has proven effective in 
the SDSS- and 2SLAQ-LRG redshift survey where a similar 
procedure were performed using the i-band photometry. Our 
star-galaxy separation algorithm only loses genuine LRGs at 
a sub-percent level and leaves «16 percent stellar contami- 
nation in the sample. 

The riz selection has been proven to work reason- 
ably well, resulting in the sample having average redshift 
2 = 0.68 ± 0.07. The redshift distribution is further con- 
firmed by the ongoing AAT-AAfi LRG project, a down- 
sized version of the VST-AAJ7 ATLAS survey, designed to 
observed several thousands LRG redshifts for photo-z cal- 
ibration and clustering evolution study. The N(z) used in 
the model prediction of ISW signal comes from » 2000 AAil 
LRG redshifts taken during a run in Jun e 2008 (Sawangwi t 
et al. 2009, in prep.) as well as those from lRoss et al.l (2008). 

The SDSS DR5 'best' imaging database contains 
800 346 photometric objects that satisfied the AAQ LRG se- 
lection rules including the necessary star-galaxy separation 
performed in the 2-band. As with the 2SLAQ LRG sample, 
objects with BRIGHT or SATURATED or BLENDED but 
not DEBLENDED flags are discarded from our sample. 
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3 THEORETICAL PREDICTION 

The secondary CMB anisotropy caused by the time-varying 
gravitional potential, $, is known as the Integrated Sachs- 
Wolfe (ISW) effect. As the CMB photons traverse such 
regions, the temparature perturbation associated with the 
time dependent potential is given by 

(1) 



A 



Ot (n) 



(n) 



To 



= -2 



dz — --— (n, z) 



dz 



where $ is the Newtonian gravitational potential at redshift 
z, n is a unit vector along a line of sight, To = 2.725 K is 
the CMB temperature at present time and zls ~ 1089 is 
the redshift at the surface of last scattering. 

The gravitational potential, $, is related to th e matter 
densi ty fluctuation via Poisson's equation (Eq. 7.14. |Peeblesl 



V 2 <E>(n, z) = 4irGa 2 p m (z) 6(n, z) 



(2) 



where a is the scale factor normalised to unity at redshift 
zero. By recalling that p C rit(0) = 3H 2 /8ivG and f2 m = 
p m (0) / pcrit (0) , the Fourier transform of the gravitational po- 
tential is 



$(fc,z) 



H \ 2 S[k,z) 



(f) 



(3) 



Unfortunately, the ISW contribution to the CMB pri- 
mary anisotropy is less than 10 percent for I > 10 and to 
make matters worse, the total anisotropy signal is dominated 
by cosmic variance at smaller I (i.e. larger an gle) where most 
of th e ISW signal is expected to be (e.g. iHu fc Scrantonl 
|2004 ). To isolate the ISW signal one must cross-correlate 
the temperature fluctuation wi th a tracer of gravitationa l 
potential projected on the sky (|Crittenden fc Turok| [l996). 
For this purpose, one can use the simple 2-point statistics 
to compute the angular cross-correlation of the temperature 
and galaxy fluctuation maps in real space, 



w g x(e) = (<5 g (ni)A T (n 2 )) 



(4) 



where n± ■ = cos#. To calculate the theoretical expec- 
tation for the real space cross-correlation, we start by com- 
puting the angular cross-power spectrum of the galaxy over- 
density and ISW temperature perturbation fields; 

-rISW 



C9T W (0 = {Sg,lm ^T,l'm' 



(•») 



Firstly, we need to expand the galaxy density fields, 
5 g (n, z), in spherical harmonics and Fourier transform them. 
For a galaxy survey with a selection function 4> g (z) and lin- 
ear bias b g (z), this is 



(2tt) 3 



dzlirji{k X )YC m (k) 
x b g (z) (j}g(z)8(k,z) 



(6) 



where ji(y) is the spherical Bessel function of the first kind 
of rank I, Yi m (k) is the spherical harmonic function and x 
is a comoving distance which is an implicit function of z 
through the relation d\ = c/H(z) dz. In obtaining Eq.|6] we 
use the or thonormality of Yi m in their expansion of a plane 
wave (e.g. IScharf et alj|l992h ; 



exp(-ifc • Ax) = 47r } i l j x (k X ) Y lm (n) Y t * m (fc) 



/ m 



(7) 



Similarly, for the ISW temperature fluctuation, by 
putting together Eq. [1] [3] and [7] this is 



d A k 
(2tt) 3 



dz ±nji{kx{z))YC m (k) 



x 3n m T 



(fee) 



) 2 - 


'8{k,z)~ 


I dz 


a(z) 



(8) 



For a flat -sky ap proximation (|Limberlll953l ). following 
lAfshordi et al.l (|2004l ) and realising that in linear perturba- 
tion theory 8(k, z) = D(z) 8(k, 0) and 

(5{ki) 8(k 2 )) = (2tt) 3 & ir ac(fei - fca) P(k) , (9) 
from Eq. El [6] and [8], Cf^ (I) can be simplified to 

4 



Cf T w (l) = 



(2i + iy 



dzP(k)W lS w(z)Wg(z) 



(10) 



Wisw{z) and W g (z) are the ISW and galaxy window func- 
tions defined as 



Wisw{z) = 3Q m T (^f) 



) 2 - 




1 dz 





and 



Wg{z) = bg{z)4>g{z)D{z) 



(11) 



(12) 



where fc m (I + l/2)/ X (z), D (z) is the linear gro wth factor 
given by the fitting formula of lCarroll et al. | l|l992i ) and P(fc) 
is the linear power spectrum at redshift zero. The survey 
selection function is given by 



<t>g(x) 



X 2 n c (x) 
f dxx 2 n c (x) 



= n (z) 



(13) 



where n c (x) is the comoving number density and n(z) is 
the normalised redshift distribution, N(z), of the galaxies 
in the survey. Finally, W„t (0) is related to the cross-power 
spectrum via the expansion in Legendre polynomials; 

oo 

21 + 1 



WgT (0) 



1=2 



4?r 



■ Pi (cos f 



)Cf T w (0- 



(14) 



We set the monopole (I — 0) and dipole (I = 1) contri- 
bution to zero, as it is done in the WMAP maps ( 32TTJ) . The 
contributions of the monopole and dipole are significan t and 
overpredict wf^ by «10 percent jCabre et alj|2006h . The 
summation in Eq. [14] converges earlier than / ~ 500 but we 
set our upper limit to I = 1000 which provides sufficiently 
stable models without sacrificing too much computing time. 
The linear power spectrum is computed using 



P(fc) = A fc" 



T 2 (fc) 



(15) 



where n s is the scalar spectral index and A is the nor- 
malisation factor with the value set by as- We use the 
transf er function, T(fc), fitting formula of lEisenstein fc Hul 
(1998). Our fiducial models assume a ACDM Universe with 
n A = 0.73, Q m = 0.27, / baryon = 0.167, a s = 0.8, h = 0.7 
and n„ = 0.95. Note that, for a fiat Universe with Q m = 1, 
the linear growth factor is equal to the scale factor, a, at all 
redshifts and Wisw(z) vanishes, hence in this case we expect 
no correlation between large-scale-structure and the CMB. 
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Figure 2. The 1° smoothed map of W-band data and galaxy number overdensity for SDSS, 2SLAQ and AAf2 LRG (Ubercal) after 
applying KQ75 and SDSS-DR5 mask. 
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4 CROSS-CORRELATION TECHNIQUE 



X {WgT,m(0]) ~ W gT (9j))] 



(18) 



Firstly, each galaxy sample is pixelis ed into equal area pixels 
on the sphere using the HEALPix (jGorski et alj|2005l ) for- 
mat, following the standard resolution and ordering scheme 
of the publicly-available WMAP5 temperature map (i.e. 
nested, res=9). The most conservative temperature mask, 
extended temperature analysis (KQ75), plus point source 
catalogue mask are then applied to the temperature maps 
(SSJ an d the pixelised galaxy distributions, discarding ap- 
proximately 30 percent of the entire sky. Additionally, in 
order to estimate fairly the galaxy background density and 
a robust cross-correlation result, the DR5 coverage mask 
including quality holes are applied to the data. We only re- 
strict the data to the most contiguous region of the NGC 
and therefore exclude the SDSS stripes 39, 42 and 43 in the 
DR5 coverage mask. After applying 'KQ75 U point source 
U DR5' mask, 516,507 out of 3,145,728 pixels (all sky) are 
admitted for the cross-correlation analysis. 

The galaxy number overdensity, S g (n), is then calcu- 
lated from the galaxy distribution maps and assigned to each 
pixel; 



S g (n) 



n g (h) 



(16) 



where n g and n g are the number of galaxies and its average 
for the sample of interest. Fig.[2]shows the W-band temper- 
ature fluctuation map and S g map for SDSS, 2SLAQ and 
AAfi LRG, smoothed with Gaussian beam of 1° Full Width 
at Half Maximum (FWHM). 

The two-point cross-correlation function at angular 
separation 9 is computed using 



w gT {9) = 



(17) 



where ft is the fraction of pixel i within the unmasked re- 
gion, hi ■ fij = cos# and At is the CMB temperature 
anisotropy measured by WMAP5 with the monopole and 
dipole contribution subtracted off. However, as we use rela- 
tively fine resolution pixels and weighting by the unmasked 
fraction does not alter our measurement, Eq. [17] is simply 
«Vr(0) = (<5 g (ni)A T (n 2 )>. 

It is a well known fact that bins in the correlation func- 
tion are correlated because the same points (or pixels in this 
case) can appear in many different pairs which are included 
in different bins, especially at large scales. To correctly esti- 
mate the statistical significance of the results, one needs to 
consider the full covariance matrix, Cy. Here, we construct 
the full covariance matrices using the jackknife re-sampling. 
In order to obtain a sufficiently stable covariance matrix, 
the jackknife subsamples of approximately twice the num- 
ber of angular bins being considered are needed. For the 
number of angular bins considered in this study, we split 
the masked temperature/overdensity map into 24 subfields 
with approximately equal area. The 24 jackknife subsamples 
are constructed from these fields, each one leaving out a dif- 
ferent subfield. The w 3 t(9) are computed for each jackknife 
subsample and the covariance matrix is 



N, 



Nj 



£ 

m — 1 



[(w gT ,m(9i) — WgT(0i)) 



where Nj = 24 in this (9i) and w g T{9i) are the 

cross-correlation measured from the mth jackknife subsam- 
ple and the average of all the subsamples in the ith bin, 
respectively. Note that the difference between w g r(9) and 
w 9 t(9) estimated using the whole sample is negligible. The 
reason for multiplying Nj — 1 is because the jackknife sub- 
samples are not independent. The statistical uncertainty for 
each individual angular bin is contained in the diagonal el- 
ements of the covariance matrix. 



5 RESULTS AND ANALYSIS 
5.1 LRG- WMAP5 

The cross-correlation results of the LRG distributions with 
the WMAP5 temperature maps using the three highest- 
frequency data plus the ILC are shown in Fig. [3] The er- 
rors are la statistical errors estimated from jackknife re- 
sampling of 24 subfields as described in 331 Generally, the 
results using different WMAP bands are in good agreement 
(within the la error) for all three LRG samples. The achro- 
matic results indicate that the contamination from effects 
such as dust, synchrotron and free- free emission which are 
frequency-dependent in nature are minimal or at least sub- 
dominant compared to our statistical uncertainties. This 
also applie s to a lesser extent to th e thermal Sunyaev- 
Zeldovich jSunvaev fc Zeldovichl [198(3) effect, although for 
the bands shown, the difference in the SZ and CMB spec- 
tral slopes is only «30 percent. However, we shall see in 
fJ7]that there is still a strong suggestion that other system- 
atic effects may still be contaminating the SDSS and 2SLAQ 
results. 

We first consider our new and higher redshift sample of 
800 000 AAQ, LRGs. This sample shows virtually no positive 
correlation with the CMB data. If anything, the data show a 
slight anticorrelation out to large scales, possibly to 9 < 1° 
(~ 30h~ Mpc at the median redshift of the sample), al- 
though the signal to noise is still low. This weak anticorrela- 
tion is observed in all WMAP5 frequency bands under study 
here (the most right column of Fig. |3J) with the exception of 
the Q band which only shows zero correlation at best with a 
possible zeropoint shift towards very large scales. As for the 
SDSS and 2SLAQ results the cross-correlation with the ILC 
map gives a systematically lower amplitude (more negative 
in AAfi case) than other bands. Given the relatively large 
scales of the null result in the AASI-WMAP5 CCF and the 
amplitude of the expected ISW signal (see Fig. H}, it would 
seem extremely unlikely that the positive correlation of the 
ISW effect could be cancelled out by the negative contri- 
bution from the thermal SZ effect. If this result is real and 
not due to some systematic effects, the implications for the 
view that the Universal expansion is accelerating, could be 
profound. 

In the case of the SDSS and 2SLAQ LRG samples, our 
results are similar to those of the previous authors who have 
analysed similar datasets. We observe marginally significant 
positive correlations in the Q, V and W bands where the 
measured u> 9 t(#)'s are similar in terms of their amplitudes 
and angular extents for each sample although the signal is 
weaker in the SDSS sample. The ILC results are slightly 
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Figure 3. The cross-correlation results of WMAP5 W , V and Q band including the ILC map (top to bottom) with the SDSS, 2SLAQ 
and A AH LRG (left to right). 



lower than the other bands in both samples but otherwise 
still within la error. Our SDSS res ults can be compared 
to the lowest redshift-bin sample of IScranton et ail (|2003l ) 
wh o used the LRGs extrac ted from the SDSS DR2 follow- 
ing lEisenstein et al.l l|200ll ) but with a much fainter mag- 
nitude limit, i < 21, and divided their samples into red- 
shift slices using photo-z. The results are similar in terms 
of amplitude but our errors are slightly smaller due to our 
larger area coverage (^7600 deg 2 as opposed to r;3400 deg 2 ) 
although their object numbers are w7 times higher than 
ours owing to the broader N(z) and fainter flux cut. The 
2SLAQ results are compar able to the 'SDSS LRG' results of 
iGiannantonio et al.l l|2008h . These authors used t he MegaZ- 
LRG photo-z catalogue of ICollister et all (|2007l ), covering 
the redshift range 0.4-0.7 with a colour-magnitude selec- 
tion similar to our 2SLAQ sample but a slightly fainter flux 
limit, idcv < 20 as opposed to 19.8. In the LRG panel of 
their Fig. 4, we see that their result has similar amplitude 



and errors (jackknife) to our 2SLAQ result. Although their 
Monte Carlo methods give somewhat larger errors than the 
jackknife estimations, the statistical significance estimated 
using errors drawn from both m ethods are very similar, 2.2- 
2.5a for their LRG catalogue. IPadmanabhan et al.l (2005) 
has also performed the analysis with a similar LRG sample 
but using the angular cross-power spectrum, Ci, making a 
direct comparison to our results difficult. The sample these 
autho rs used is somewhat similar to the lEisenstein et al.l 
(2001) selection but with the flux cut as faint as 2SLAQ 
in 'CutlP, resulting in a redshift distribution similar to our 
SDSS and 2SLAQ LRG samples combined, although they 
limited the redshift of the sample to 0.2 < z < 0.6 us- 
ing their template-fitting ph oto-z. The positive correla tion 
is detected at 2.5<r, similar to IGiannantonio et al.l (|2008l ) al- 
though the sample they used only covers half as much sky. 
We conclude that our analyses are broadly reproducing pre- 
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Figure 4. The LRG- WMAP5 cross-correlation results using W-band and ILC map compared to the theoretical predictions (red solid 
lines), assuming the standard ACDM and the galaxy linear bias (b g ) of 2.10, 1.99, 2.2 and 2.1 for SDSS, 2SLAQ, AAC-LRG and the 
combined sample respectively. The stellar contamination correction for each sample has been applied to the corresponding model. In the 
'Combined' panels, the cross— correlation results of the quadrature— error weighted mean of the three LRG samples are compared to the 
average model predictions. 



vious results in these 0.25 < z < 0.6 LRG redshift ranges, 
both in terms of their amplitude and statistical significance. 

5.2 Comparison to Models 

Fig. [4] shows the comparison of our results to the theoret- 
ical expectation as described in ij3] The galaxy selection 
functions used in construction of these models are given by 
the normalised N(z) of the sample as shown in Fig. Q] (see 
also i|2.2[) . The galaxy bias in the model is estimated from 
the angular autocorrelation function, w g (8,z), of each LRG 
sample relative to the underlying dark matter clustering, 
b 2 g (z) — £ g (r,z)/t; m (r,z), where we assume the linear scale- 
independent bias and measure its value at large scales (~10 
h~ 1 Mpc). The matter { (r, z) is estimated for the same fidu- 
cial cosmology as described in §5\ and then projected onto 
the sky using the corresponding n g (z). This gives an un- 
biased prediction which can be compared to the measured 
w g (8, z) and allows b g (z) to be extracted from their relative 
amplitudes (see Sawangwit et al. in prep for the full detailed 
analyses). Note that, we assume non-evolving bias and de- 
note the bias estimated from each sample as the bias at the 
corresponding average redshift which is reasonable, given the 
narrow redshift ranges of our samples. The galaxy bias mea- 
sured in this way can also be viewed as an effective value 
for each sample. The models shown in Fig. [4] use b g (z) of 
2.10±0.04, 1.99±0.02 and 2.20±0.02 for the SDSS, 2SLAQ 
and AAf^samples, respectively. These values are taken from 
Sawangwit et al. (in prep) and are c ompatible with t he va l- 
ues measured by other aut hors, e.g. iTegmark et al] |2006), 
iPadmanabhan et al.1 (|2007l ) whose 6 9 (0.35) = 1.9 ±0.07 and 
6 9 (0.55) = 1.85 ±0.05 as compared to our SDSS and 2SLAQ 



LRG, respectively. The b g value of lTegmark et all (2006) was 
measured from a sample of z ~ 0.35 LRGs similar to what 
we call SDSS LRG sample here but without the bright limit 
cut (see i|2.2.1fl hence allows under-luminous objects and 
main galaxies into their sample. And as a result we expect 
their bias to be somewhat lower than ours. 

As emphasised earlier, the AAQ LRG sample shows 
no positive correlation with the WMAP5 data and perhaps 
even a slight negative correlation. We then combined the W- 
band data between 12'-120', and found the amplitude of the 
CCF and its jackknife error (la) is -0.07 ± 0.2 /xK. This is 
consistent with the null hypothesis (only ~ 0.4cr deviation) 
and rejects the ISW signal expected in the standard models 
at pb 1.9a or at 5% significance after the stellar contamina- 
tion has been taken into account in the predicted signal (see 
N6.3p . Performing a similar statistical analysis on the cross- 
correlation results using the ILC map gives a slightly higher 
significance of rejecting the standard model ISW hypothesis 
(2.2a, see Table 0). 

Additionally, to test that the zero correlation in the 
AAQ results is not due to its faint limit making the sam- 
ple incomplete, we have cut the faint limit of the sample 
back in steps of 0.25 mag to 20.0 (see ij6.2l and Sawangwit et 
al. 2009). The amplitude of the CCF between 12'-120' for 
i < 20.25 (denoted by AAfi* in Table© is -0.1 ±0.2 for W- 
band data and -0.2±0.21 for the ILC map. The ISW model 
prediction is then re-computed taking into account the cor- 
responding n(z) and linear bias, including the correction for 
stellar contamination at the same level as the main AAQ 
sample. The significance of rejection of the standard model 
for the i < 20.25 AAQ sample is slightly higher than that 
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Table 2. The significance tests of the cross-correlation results using the WMAP W-band data and ILC maps. The measurements are 
tested against the expected ISW prediction in the standard ACDM model and null result hypothesis. Column 5 gives the amplitudes 
and la jackknife errors of the data binning between 12'-120'. Column 6 gives the significance of the deviation of the value in column 5 
relative to ISW/null signal hypothesis. 





Sample 


z 


Number 


ba(z) 


w gT (12' - 120') 


Deviation significance 












fiK 


(ISW.mill) 


W-band: 


SDSS 


0.35 


106,699 


2.10 ± 0.04 


0.25 ± 0.33 


(l.Ocr, 0.8<r) 




2SLAQ 


0.55 


655,775 


1.99 ±0.02 


0.34 ± 0.21 


(0.2ct, 1.6a) 




AAf2 


0.68 


800,346 


2.20 ±0.02 


-0.07 ±0.20 


(1.9cr, 0.4(7) 




AAQ* 


0.67 


375,056 


2.37 ±0.03 


-0.10 ±0.20 


(2.2<7, 0.5a) 




Combined 


0.60 


1,562,820 


2.10 ±0.01 


0.15 ± 0.17 


(1.0(7,0.9(7) 




Weighted mean 








0.14 ±0.14 


(1.3(7, 1.0(7) 


ILC map: 


SDSS 
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2.10 ±0.04 


0.19 ± 0.33 


(1.2(7,0.6(7) 
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1.99 ±0.02 


0.27 ±0.22 


(0.5cr, 1.2cr) 




AAf2 
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-0.18 ±0.22 
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(2.5(7, l.Ocr) 
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0.60 
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2.10 ±0.01 


0.07 ±0.17 


(1.4(7,0.4(7) 
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Figure 5. Top: The W-band cross— correlation results of the com- 
bined sample (solid circles) compared to the quadrature-error 
weighted mean of the three LRG sample (diamonds). Also shown 
are the standard model predictions by taking a weighted mean 
(solid line) of the models of three LRG samples and for the com- 
bined sample (dot-dash line). Bottom: Same as above but for the 
ILC map rather than W-band data 



of the main AAil sample, at 2.2a and 2.5a for IF-band and 
ILC map, respectively. 

The measured w g r for the SDSS LRG agrees reason- 
ably well with the theoretical expectation at angles < 30' 
although not at high statistical significance. However, the 
same cannot be said for the angle beyond this scale and up 
to ~ 600' where the cross-correlation appears to be less than 
the expected signal although still not at high significance. 
One may be inclined to conjecture that this could be due to 
the negative contribution coming from the thermal SZ effect, 
but at this redshift 100' corresponds to ~ 20 /i~ 1 Mpc which 
would be too large a scale to be caused by hot gas in galaxy 
clusters. Although the clusters do cluster among themselves, 
the contribution to any extend SZ effect is likely to be small 
<|Mvers et al.ll2004l ). Besides, there is no physical reason why 
should SZ effect only affect the highest redshift sample. The 
most likely explanation for this appears to be a statistical 
fluctuation which means that our SDSS LRG measurement 
rejects neither the ISW expectation nor the zero correlation 
at more than « lcr significance level. If we bin the data 
in the angular range 12'-120' into a single bin, the corre- 
lation amplitude and its jackknife error (la) is 0.25 ± 0.33 
/xK which deviates from the null result hypothesis by only 
0.8a and from the standard model by l.Ocr. For the 2SLAQ 
case, as in other studies, the positive cross-correlation sig- 
nal agrees very well with the expected ISW signal in the 
standard cosmology in terms of its amplitude and angular 
extent. Nevertheless, the 2SLAQ sample's rejection of the 
null result is still only at the 1.2-1.6<r significance level (see 
Table 



5.3 The Combined LRG sample 

We shall now consider the cross-correlation of the combined 
LRG sample with the CMB data. In our first method of 
combining the three LRG samples we shall treat these as 
three independent surveys and then test this assumption by 
presenting the cross-correlation result for the combined 1.5 
million LRG sample, complete with its own direct jackknife 
error analysis, to check that they agree. 
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First, the three CCF's of the SDSS, 2SLAQ and AAfi 
samples are combined by weighting inversely in quadrature 
according to the statistical errors of each sample (see bottom 
right panels of Fig. [4] and also Fig. [5j . The model (red solid 
line in Fig. [4| is estimated by taking an average of the ISW 
models of the three LRG samples. We find that the rejection 
significance is 1.3a for the standard ISW model and I.Oct 
for the null result in the W band. In the ILC band the 
significance of the rejection of the ISW model rises to 2. Oct 
and the significance of rejection of the null result reduces 
to 0.5cr. Table [2] gives the summary of all the significance 
tests performed. We conclude that while the ISW standard 
model is still consistent with the CCF result from the three 
combined, weighted LRG samples it is now more consistent 
with the null result due to the inclusion of the AAQ data. 

Second, for comparison, we also present the results 
of cross-correlating the combined LRG sample with the 
WMAP5 data i.e. we now treat the combined sample of w 
1.5 million LRGs as a single sample for cross-correlating 
with, in turn, the WMAP5 W and ILC CMB data. A full 
jackknife error analysis was carried out in the same way 
as for the individual samples. We expect the results to be 
similar to the weighted combination of the three samples' 
CCF's as presented above. Fig. [S] shows the comparison be- 
tween these results. The models for the combined samples 
are computed following the procedure described in Sj3] ass- 
suming the linear galaxy bias (given in Table 0) estimated 
from the angular autocorrelation function and N(z) of 1.5 
million LRGs. Table [2] again shows the significances of re- 
jection of the standard model and the null results. We see 
that the observational results in both cases are very similar. 
For both bands, the significances are given in Table [2] The 
results are again very similar to those where the weighted 
mean was adopted. The cross-correlation results are again 
as consistent with the zero correlation as they are with the 
standard ISW model for the W band. The ILC band again 
more signifcantly rejects the ISW model than the null result. 

Clearly the preference for the null result over the stan- 
dard model prediction depends on the accuracy of the new 
AAQ. result. We test the robustness of the AA57 result in 35] 

5.4 x 2 fits 

For completeness, we also quantified the goodness-of-fit of 
our measurements to the expected ISW signal or null result 
hypothesis by calculating the chi-square, \ 2 , which uses the 
normal size bin as shown in Fig. [4] and takes into account 
the possible correlation of the bins through the use of the 
covariance matrix (©. The x 2 is given by 

* = E ^W(*i)-«& W (*i)] • Kt(^)-< t w (^)](19) 

where Cr. is the inverse of covariance matrix, w g T{di) is the 
measured angular cross-correlation and is the theoret- 

ical expectation assuming the standard ACDM cosmology 
(see Fig. 2} which can be replaced by zero when testing the 
zero correlation hypothesis. Using the galaxy linear bias, b g , 
and N(z) for each sample as as mentioned in EI5.2I the \ 2 
tests were performed for the angular bins between 12'-120', 
inclusively. The lower limit is set approximately to the best 
WMAP5 resolution in the W band (» 12'). 



The significances obtained from the \ 2 method gen- 
erally confirmed the results using the 12'-120' bin, espe- 
cially those of the main LRG samples. For example, assum- 
ing standard model parameters, the SDSS-W band results 
give x 2 = 19.4 for the predicted ISW signal and 17.7 for the 
zero correlation hypothesis. For the 2SLAQ results, using 
the standard model gives \ 2 ~ 13.2 and relative to the 
null result gives \ 2 — 11.5. These \ 2 were computed for 
6 degrees-of- freedom (d.o.f). Using the x 2 distribution, the 
SDSS results deviate from the ISW model and null result 
at 4 and 7 percent statistical significance, respectively. The 
2SLAQ results agree with the ISW model with the reduced 
chi-square, x 2 /d-o.f., of order of unity and reject the zero 
correlation hypothesis at 1.5er significance. The AAQ re- 
sults gave x 2 = 11-7 and x 2 = 4.4 for the ISW model and 
null correlation respectively, corresponding to the chances of 
agreement of 7% and 62%. These all agree reasonably well 
with the large-bin significances presented in Table [2] How- 
ever, the similar x 2 significance tests of the combined sample 
and some ILC individual samples did not perform very con- 
sistently, occasionally giving pathological results and poor 
agreement with the 12'-120' bin and this is why we have 
only quoted the simpler, single large-bin significances in Ta- 
bled 



6 ROBUSTNESS TESTS 

Given that the AAQ LRGs comprise a new sample, there is 
no previous measurement that can be directly compared to 
our own. We now present the result of tests we have done in 
order to check the robustness of our new result. 



6.1 Random realisations and simulated CMB 
Maps 

Here we generate 100 random realisations for each of the 
sample. Each realisation has the same number density as 
the sample it tries to mimic. Note that these random re- 
alisations are unclustered. The results are shown in Fig. [5] 
The jack-knife errors that we use are seen to be much larger 
than the standard deviation of the random catalogues (inner 
green shaded region). This is expected because the random 
catalogues are unclustered unlike the LRGs. The means of 
these random realisations are consistent with zero and show 
no sign of bias except perhaps at the smallest scales of the 
SDSS sample. 

We have also made simulated CMB temperature 
anisotropy maps and cross-correlated these with the three 
LRG samples. A simulated CMB map is generated as a re- 
alisation of random gaussian fields on a sphere with the 
fluctuation characterised by WMAP5 best-fit power spec- 
trum. The simulated maps are also convolved with a Gaus- 
sian beam with FWHM similar to the WMAP W-band, i.e. 
12'. 6. The cross-correlation results are shown in Fig. [6] The 
standard deviation of 100 CMB random realisations (outer 
grey shaded region) are roughly consistent with our jack- 
knife estimates especially at small and intermediate scales 
but somewhat larger at large scales. 
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Figure 6. The cross-correlation results (diamonds) of three LRG samples and their jackknife errors (lcr) compared to the results of 
using 100 random realisations of each LRG sample (inner green shaded region) and 100 simulated CMB maps (outer grey shaded region). 
The shaded area signify a standard deviation in the measurement of 100 realisation for each case. Note that the means of these random 
realisations are consistent with zero as can be seen from their symmetry about the zero line. The solid line is again the theoretical 
prediction of the ISW signal in standard ACDM. 
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Figure 7. The cross-correlation of the AAQ LRG to W-band 
data using the original SDSS photometry (diamond) compared 
to the measurements using 'ubercalibration' (dot-dash line), the 
stripe weighted (dotted line) and when the data is restricted to 
the region where galactic extinction in i-band less than 0.1 (dot- 
dot-dot-dash line) and 0.08 (long dashed line) magnitude. 



6.2 Photometry test 

Next, we look to see if the AAf2 cross-correlation measure- 
ment is robust by comparing the resu l t from the SDSS 'uber- 
calibration' of iPadmanabhan et al. (2008) with that from 
the standard SDSS calibration. Fig. [7] shows that the re- 
sults are stable to whichever calibration we used. We fur- 
ther looked for systematic effects in the original photometry 
by weighting SDSS stripes according to their overall num- 
ber density. The physical motivation for this arises from the 
SDSS observing strategy and the fact that a slightly differ- 
ent calibration for different nights could affect the source 
density as a function of the SDSS stripe, given our faint 
limit. We observe a hint of these variations although not at 



a high level and use these to correct the source densities in 
each stripe as mentioned. However, such variations seems 
to be weaker when using the ubercalibration as opposed to 
the standard one. The result of weighting according to the 
stripe number density is shown in Fig. [7J and again the re- 
sult appears robust when this filter is applied to the original 
data. 



Although we work at a relatively high galactic latitude, 
it is possible that in some regions of the sky, high galactic 
dust obscuration could result in lower detections of faint ob- 
jects. Furthermore, that same dust obscuration patch could 
be a source of contamination in the CMB data in the sense 
that the temperature in that particular region could be sys- 
tematically raised by the dust emission and hence results in 
a false anticorrelation. To test this, we exclude the region 
where the extinction is greater than 0.1 mag in the i-band 
which discards ~ 15 percent of the data. We observe no 
difference to our main results, even when a more aggressive 
limit, ^extinction < 0.08 (23 percent discard), is applied (see 
Fig. 0. Note that when similar tests are performed using 
extinction in the SDSS r-band instead, we again obtain re- 
sults which are consistent with those presented in ij5]for all 
three LRG samples. 



We then cut back the i-band limit of the AAS1 sample 
in 0.25 mag steps from i = 20.5 to i = 20.0 while keeping the 
other conditions the same. These results are compared with 
the result at i < 20.5 in Fig. [8] Again the results appear ro- 
bust. We have also made tests of the single epoch SDSS pho- 
tometry using deeper Stripe 82 (lAbazaiian e t al. 200 °j) and 
the W illiam Herschel Deep Field (WHDF, iMetcalfe et~ai1 
l200ll ) data. Both these comparisons showed that the SDSS 
photometry in r, i and z bands showed good agreement with 
the deeper data until the errors showed a significant increase 
beyond the limits r = 22.0, i = 21.0 and z = 20.2. 
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Figure 8. The cross-correlation of the AAQ LRG to W-band 
data compared to the measurements using the same colour-colour 
selection sample but with brighter faint-limit cut, i^cV < 20.25 
and «deV < 20.0. Only the theoretical expectation of the full 
(solid line) and i^cV < 20.25 (dash-dot line) sample are shown. 
The data points are shifted slightly for displaying purposes. 
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Figure 9. The object numbers per pixel as a function of galactic 
latitude, b. Recall that we use equal are a 49 arcmin 2 ) pixels 
with res= 9 resolution scheme fHEALPix. lGorski et al.l2005h . The 
SDSS number has been multiplied by 10 to extend the plot range. 



6.3 Star— galaxy separation 

We noted in fj2]that the stellar contamination in our AASl- 
LRG sample could be as high as 16 percent. Care should 
therefore be taken when analysing this dataset. We obtained 
this contamination fraction using the inf ormation learned 
from the AAil-LRG spectroscopic survey IjRoss et al.ll2008i. 
Sawangwit et al. in prep), by imposing a star-galaxy sep- 
aration in the z-band similar to the method applied in the 
SDSS- and 2SLAQ-LRG redshift surveys using the i-band. 
Our optimised star-galaxy separation procedure selects ob- 
jects with Zp S f — Zmodoi > 0.53 + 0.53(19.0 — z mo dei) which 
only loses genuine LRGs at a sub-percent level and leaves 
~16 percent stellar contamination in our sample, as quoted 
eariler. 
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Figure 10. The cross-correlation of the W-band data and the 
2SLAQ LRG when a sample of random realisation of £Sl6 per cent 
is added to the LRG catalogue (diamonds). The results using 
original 2SLAQ sample and when multiplied by 1 — / s are shown 
as the dot-dashed and dotted lines, respectively. The long-dashed 
line shows the result when the 16 percent added contaminants 
are replcaed by red stars. The result of cross-correlating the ILC 
map with the 16 percent-red star contaminated 2SLAQ sample 
is also shown (solid line). 



The effect of stellar contamination distributed at ran- 
dom in the sample is simply a dilution of the over/under 
density hence reducing the autocorrelation amplitude of the 
sample by (1 — f) 2 and the cross-correlation by a factor of 
(1 — /) where / is the fraction of the contamination. This 
is particularly true if the contamination is distributed uni- 
formly at random in the sample. However, if there is some 
spatially dependent variation of the number density, a fur- 
ther systematic effect could arise through this process. To 
test this, we first check to see if there is a trend of the number 
density as a function of galactic latitude as one might ex- 
pect for stellar contamination. Although a slight such trend 
is observed, it is at no more than the levels observed in 
the SDSS and 2SLAQ samples (see Fig. [9} whose stellar 
contamination fractions are approximately 1 and 5 percent, 
respectively. Next we restrict the data to the high galactic 
latitude regions, namely b > 40° , 50° and 60° . The results 
are in good agreement with our main results for all three 
LRG samples up to b > 60° where the cross-correlations 
become noisy due to the 75 per cent reduction in the sample 
sizes. 

To simulate the effect of the stellar contamination on 
the LRG-CMB cross-correlation, we have introduced a set 
of random realisations into the 2SLAQ sample. The result 
is presented (diamonds) in Fig. [10] along with the cross- 
correlation of the original 2SLAQ sample (dot-dashed line) 
and the result of reducing its amplitude by a factor of 
(1 — 0.16) (dotted line). Furthermore we would like to check 
for any effects that may arise from possible large scale clus- 
tering of the stars. This is done by adding a sample of red 
stars to the 2SLAQ LRG sample at the 16 percent level. 
The stars are selected with similar colour-magnitude crite- 
ria to that of the AAS1 LRGs and should mimic the angular 
distribution and properties of the stellar contaminants seen 
in the sample. The result is shown in Fig. [10] (long-dashed 
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Figure 11. Left: The AAH LRG-WMAP5 cross— correlation of the 9% stellar contaminated sample (asterisks) compared to the main 
AAf! sample used in our study (diamonds). Right: The corrected autocorrelation functions of the 9%— and 16%-contaminatcd samples 
(asterisks and diamonds). These are compared to the results of limiting the 9%-contaminated sample to the regions with galactic latitude 
higher than > 40°, 50° and 60°. The dashed-li ne and shaded re gion is the w gg {9) and lcr error inferred from the w p (cr) measured from 
~ 400 spectroscopically confirmed AAfl LRGs llRoss et al.ll2008f) . 



line). This test should also reveal any possible effects on the 
io 9 t's due to (if any) correlation between these stars and 
the CMB. We found the 16 percent-red star contaminated 
2SLAQ result to be consistent with the dilution of randomly 
distributed contaminants case. The result is also consistent 
with the cross-correlation with the foreground reduced ILC 
map (solid line), further confirming that our result is not 
affected by any star-CMB cross-correlation. Note that the 
significance test presented in Table [2] has already taken into 
account such an effect by multiplying the ISW model by a 
factor of (1 — 0.16). The significance of the AAfi sample's 
rejection of the standard model ISW prediction is therefore 
robust against the stellar contamination discussed here. 

We next attempt to reduce the stellar contamination 
fraction by imposing a more aggressive star-galaxy separa- 
tion cut which result in nearly halving the number of genuine 
AAfi LRGs. The cut is a combination of the fitted 'de Vau- 
couleurs' radius as a function of ZdeV magnitude and the cor- 
relation between the 'de Vaucouleurs' and fiber magnitudes 
in z-band. This reduces the contamination to ~9 percent. 
Fig. [TT] (left panel) shows the cross-correlation of this sam- 
ple with the W-band data which is in good agreement with 
our main result. 

The contamination fractions of these samples are veri- 
fied by their angular autocorrelation functions, w gg (8). The 
corrected w gg (8) is shown in right panel of Fig. [11] This 
is again in good agreement with the 16 perce nt contami- 
nated sample and consistent within ~ lcr of the Ros s et all 
(2008) power-law fit to the semi-projected correlation func- 
tion, w p (a). Note that we only expect the agreement in the 
range r ~1— 15 h~ 1 Mpc where a single power-law is a good 
fit to the data. The measured w gg 's are also consistent with 
the results when restricting galactic latitude to greater than 
40° , 50° and 60° . We believe the slight discrepancy with the 
w p (a) is due to the noisy measurement from the small num- 



ber of spectroscopically confirmed LRGs used in iRoss et"al1 
(2008) and not caused by the under-estimation of the con- 
tamination level as demonstrated by our two independent 
approaches for star-galaxy separation. 

Even if the contamination fraction is under-estimated, 
the effect of an increased (uniform) stellar contamination 
would be to increase the ISW model amplitude when the bias 
value from the LRG autocorrelation is corrected upwards to 
obtain the true bias value. This upwards shift in the ISW 
model would then be exactly cancelled by the downwards 
correction to account for the dilution of the cross-correlation 
signal due to stellar contamination. 

We conclude that despite the faint magnitude limit, and 
moderate level of stellar contamination (~ 16 percent) our 
ISW results for the AAS1 LRGs seem robust to the tests we 
have made and the SDSS data seem accurate enough to sup- 
port this ISW analysis. Up to this point, we have therefore 
found no explanation in terms of a systematic effect for the 
low AAQ- WMAP5 cross-correlation result. Next, we shall 
perform a similar analysis on some of the large-scale tracers 
whose ISW effect has been previously claimed in order to 
test our methodology and look for other possible systemat- 
ics in these samples. 

6.4 SDSS galaxy- WMAP5 

We next cross-correlate galaxies extracted from SDSS DR5 
using r-band magnitude limits. The objects are photomet- 
rically classified as galaxies by the SDSS reduction pipeline. 
We subsample the galaxies in three magnitude ranges, 
namely, 18 < r < 19, 19 < r < 20 and 20 < r < 21, 
where all the magnitudes are galactic extinction corrected 
model magnitudes. The subsamples contain approximately 
2, 6 and 16 million s objects, respectively. This is the same as 
ICabre et all (|2006l . C06 hereafter) but covering w 20 per cent 



© 2009 RAS, MNRAS 00Q.[lTfl9l 




Figure 12. The cross— correlation of W-band data with the r-band selec ted galaxies. The samp le magnitude ranges are as indicated in 
the plots. The ISW model prediction is shown for each sample, assuming iDodelson et aD (12002^ 1 n(z) with the measured bias of 1.2, 1.1 
and 1.2 for 18 < r < 19, 19 < r < 20 and 20 < r < 21 sample, respectively. 



more area and we use WMAP5 rather than WMAP3. A 
similar r-band sele c ted galaxy sample was also used by 
iGiannantonio et alJ (|2008l . G08 hereafter) although these 
authors use 'ubercalibration' photometry rather than the 
original one and limit the sample photo-z to redshift between 
0.1-0.9. The ISW effect has been claimed to be detected in 
these samples at moderate significance level by both C06 
and G08 although their results do not agree with the former 
having twice as much positive cross-correlation between the 
CMB and the r-band selected galaxy sample. As a result 
C06 need to fit their result with higher Qa , for a galaxy bias 
b = 1.0. 

For the cross-correlation analysis, we proceed in the 
same manner as with the LRG samples. To compute 
the ISW mode l , we use the n(z) distributions following 
IDodelson et"aH i|2002D . The average redshifts inferred from 
the n(z) are estimated to be approximately 0.17, 0.24 and 
0.33. We then follow our procedure for the LRGs and obtain 
the galaxy linear bias from the measured amplitude of the 
galaxy 2-point autocorrelation function for each subsample. 
We obtain the values b g — 1.2, 1.1 and 1.2 for the sample 
with r-band magnitude limit of 18-19, 19-20 and 20-21, re- 
spectively, in agreement with the measurements of C06 and 
G08 whose b g ml— 1.2. The cross-correlation measurements 
and the theoretical models are presented in Fig. 1121 

We marginally detected the correlation between the 
CMB data and all the r-band selected subsamples. We shall 
now compare the 20 < r < 21 result in Fig. ll2l to Fig. 2 (top) 
of C06. Our result is lower by a factor of ~ 2 but very close 
to the re-analysis of the SDSS r-band data of G08 who also 
found a factor of two discrepancy with C06. After their dis- 
cussions, the two groups found that the discrepancy is due to 
an extra quality-cut imposed on the data by C06, namely, 
r-band magnitude error less than 0.2 mag. We regard the 
factor of two rise in the amplitude of the cross-correlation 
after this small change in the magnitude error limit simply 
as symptomatic of the statistical fragility of the result. We 
conclude that our re-analyses of these data agree well with 
the standard ACDM predictions although the significance of 
rejection of the null result is still only «l-2a. 



6.5 NVSS- WMAP5 cross-correlation 

To test our methodology further, we performed a cross- 
correlation analysis of WMAP5 with radio sources fr om the 
NRAO VLA Sky Survey (NVSS; ICondon et al.|[l99Sh which 
has been previously used by various groups for ISW stud- 
ies. The NVSS sample comprises about 1.8 million radio 
sources detected to a flux limit of » 2.5mJy at 1.4GHz. 
The NVSS covers the entire sky higher than —40° decli- 
nation (~80 percent of the sky). Intere stingly, the previ- 
ous study of iBoughn fc Crittenden! l|2002l ) found no correla- 
tion of these sources with the Cosmic B ackground Explorer 
(COBE) CMB map but a later study bv lNolta et"al] (|2004[ ) 
found a positive correlation with the first-year WMAP data 
which they claimed to be the evidence for > at 
95% confidence, assuming a flat CDM cosmology. The re- 
analysi s of the NVSS-CM B correlation by G08 also con- 
firmed iNolt a et alJ (|2004l ) results at approximately the 3<r 
significance level. 

For the cross-correlation analysis we restrict the data 
to the declination, 8 ^ —37° where the survey is the most 
complete. We then applied the masking and pixelisation pro- 
cedure described in §3] but for this sample we shall use 
lower resolut i on (r es=6 as opposed to res=9) HEALPix 
iGorski et all (2005) scheme to reduce the computing time 
because of the much larger sky coverage of the NVSS. We 
checked that the measurements using different resolutions 
do give the same results in terms of amplitudes and statisti- 
cal uncertainties (§5J. The higher resolution (res=9) result 
shall be discussed in this section but for the purpose of the 
systematics test in iJT] we shall present the results using 
res=6. 

Boug hrTfe Crittenden! l|2002l 'l noticed a number density 
trend with the declination w hich affected their autocorrela- 
tion measurement. Following iNolta et al.l (|2004T ). we applied 
a correction for this by splitting the sample into sin 8 strips 
of width « 0.1 and scaling the galaxy numbers in pixels be- 
longing to a particular strip by the ratio of global mean to 
the strip mean. The cross-correlation procedure is then car- 
ried out as outlined in i|4]but the statistical uncertainties and 
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Figure 13. The cross-correlation of the NVSS sources to the W- 
band data . The ISW prediction (red solid line), assumes linear 
bias of 1.5 llBoughn fc Crittendenl2002l;lGiannantonio et al.l2008t) 
and n(z) derived from IPvmtorT'fcPelilcockl (| 1 9901 1 radio source lu- 
minosity function (mean-z model 1). 

covariance matrix are now estimated from approximately 20 
equal-area jackknifes rather than 24. The result using W- 
band data along with the standard model ISW prediction 
(red solid line) is presented in Fig. [13] 

The ISW predictions for the NVSS sources are com- 
puted using the number-redshift distribution, n(z), derived 
from the radio source luminos ity function (mean- z model 
1) of lDunlop fe Peacock! l|l990Ti . The median redshift esti- 
mated from such n(z) is ~ 0.8 with a tail extending out to 
z ~ 3. We assume the source bias, b, of 1.5 as measured 
by a number of authors (e.g. iBoughn fe Crittenden! [20021; 
iGiannantonio et alJliooih . 

Fig. [T3] shows that we find a marginally positive cor- 
relation similar to the prediction of the standard model 
at scales smaller than « 5° at ~ 2a significance. Our re- 
sult can be directly compared with that of G08 who, like 
us, cross-correlate the source number fluctuations to At 
as opposed to th e source number per pixel approach of 
iNolta et al.l l|2004h . Similary, G08 observed a good agree- 
ment between their measurement and the standard ACDM 
model which also starts to breakdown at ~ 5°. We take 
this agreement as a further indication of the robustness of 
our cross-correlation methodology and analyses. In the next 
section, we shall further test the NVSS-WMAP5 result for 
contamination by systematic effects. 



7 CMB SKY ROTATION TEST 

Here we shall perform an ad ditional test f or sy stem- 
atics, similar to that used by I Myers et al.l l|2004h and 
iBielbv fe Shanks! (2007) for testing their detection of the 
SZ effect, particularly in checking the reality of a large scale 
temperature decrement around galaxy groups and clusters. 
We follow these authors and rotate the WMAP maps around 
the galactic pole in the clockwise direction, each time adding 
40° to galactic I. There is an area very close to the pole where 
there is less movement from the rotation, but given that we 
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Figure 14. The cross-correlation of the three LRG samples (top), 
r-band selected galaxies (middle) and the NVSS sources (botom) 
to the rotated W-band data in our rotation test (see text for 
more details). Note that for the top two plots, the points have 
been shifted slightly in the x-axis for clarity. 



use a 40° shift the effect of this slight non- independence is 
small. We have checked that if we cut out the circum-polar 
region down to galactic latitude b = 75° our results are un- 
affected. 

The CMB masks (KQ75 plus point source) are rotated 
with the temperature maps to ensure that the contaminated 
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regions are excluded from both galaxy and temperature fluc- 
tuation maps. The SDSS DR5 mask is then applied to the 
data in the case of LRG and r-band selected samples. The 
cross-correlation is performed using the W^-band data fol- 
lowing the procedure described in Sj4] We use the cross- 
correlation results between 12' < 8 < 120' where the dif- 
ference between the ISW and null result is at its maximum 
as in SJU The cross-correlations are then performed at eight 
40° intervals. 

7.1 LRGs 

The cross-correlation measurements are presented in Fig. [T3] 
(top panel). The errors shown are jackknife errors (la) and 
as expected they are similar at all rotation angles which 
makes the data points straightforward to compare. For the 
SDSS sample at z=0.35, there is a higher positive point at 
rotation angle 40 degrees. For the 2SLAQ sample at z=0.55 
the points at rotation angles 160 and 240 degrees are more 
negative than the zero degree point is positive. There is no 
reason to expect anything other than a null result at any ro- 
tation angle other than zero. Therefore, based on this rota- 
tion test the significances are now reduced to the « 12 — 25% 
level, suggesting that systematics as well as statistical errors 
may be affecting the data. 

7.2 SDSS galaxies 

We also applied the same test to the ISW results using three 
SDSS r-band selected galaxy samples of 18 < r < 19, 19 < 
r < 20 and 20 < r < 21. The results are shown in Fig. [14] 
(middle). Again we see that there are rotation angles that 
show more significant non-zero detections than at the zero 
degree rotation angle. We see that at 40° rotation angle, the 
results are very negative in all three samples. At the rotation 
angle of 200°, the results are more positive than the zero 
degree rotation, again in all three samples. As for the LRG 
samples, this means that the significance is reduced to a 
marginal ~ 10% level and the results suggest that systematic 
effects as well as statistical errors may be contributing to the 
apparent ISW detection at zero degrees rotation angle. 

7.3 NVSS radio sources 

We then applied the same test to the NVSS- WMAP5 cross- 
correlation result, (see Fig. 1141 bottom). This time the point 
at rotation angle 280 degrees is more positive than the point 
at zero degrees. As with other samples, the jackknife errors 
on all the points are similar so this comparison is fair. Again 
we conclude that systematic effects may be contributing to 
the apparent ISW detection which explains the reduction in 
statistical significance to > 10% from the rotation test. 



8 DISCUSSION 

Given the consistency of the AAfl and the combined LRG 
results with the zero correlation, we now discuss whether 
there is any contradiction between our conclusions and those 
of other authors. In particular, we discuss the results of G08 
who claim a 4.5a ISW detection from the combined analyses 
of several large-scale tracers. These tracers include some of 



the LRG samples. They also include NVSS radio sources. 
The most significant detection in their Table 1 is from the 
NVSS at 3.3<r. Their LRG analysis gives 2.2cr for a sam- 
ple roughly equivalent to our 2SLAQ LRG sample. These 
compare to 1.6a for our 2SLAQ samples. For the NVSS we 
find a 1.8<r result. Their SDSS galaxy sample gives 2.2cr 
equivalent to our combined SDSS i — band limited sample 
which gives ~ 1.3ct significance. Thus our significances ap- 
pear lower than those of G08, particularly for NVSS. This 
discrepancy increases when we consider the rotation test. In 
the rotation test of the NVSS sample, 1 out of 8 points has 
higher amplitude than zero rotation measurement which is 
only w 1.5a significance. For the 2SLAQ case, this gives 1-2 
out of 8 points which is equivalent to 1.2-1. 5a significance. 
The SDSS galaxy gives 2 higher (or lower) points in 8 or 
» 1.2a. 

However, G08 also have differences between their two 
methods of assessing the significance of their results. Their 
Table 1 assumes the hypothesis of the standard ACDM 
model to obtain a maximum likelihood amplitude, A, and an 
associated error from their data. This error is different from 
the error that can be inferred from the x 2 statistic in their 
Table 2 which tests the null result hypothesis. For example, 
their LRG result is 0.4er significant from Table 2 whereas it is 
2.2a significant from Table 1. Their SDSS galaxy sample re- 
jects the null result hypothesis at 1.3<r significance from the 
X 2 statistic, again lower than their Table 1 at 2.2cr. Also the 
NVSS only reject the null result at 1.3a rather than 3.3<r. We 
assume that these differences may be due partly to different 
null hypotheses and partly due to different methodologies. 
Certainly, the levels of significance in their Table 2 are lower 
and more in line with what our rotation tests show, i.e. 1-2 
higher (or lower) points in 8 or 1.2-1.5<r. It remains to be 
seen for the other samples in their Table 1 and 2, if the same 
pattern applies with the maximum likelihood significances 
in Table 1 being higher than the x 2 significances in Table 2. 
We conclude that the rejection of the null result by their x 2 
test may be more consistent with what we have found than 
the results in their Table 1. Indeed their x 2 summed from 
all surveys is 67 on 74 degrees of freedom which is hardly a 
significant rejection of the null result and can be compared 
to our overall rejection of the null result in our Table 2 of 0.5 
to la. Therefore as long as we refer to the x 2 test of G08, 
there seems to be no inconsistency with our estimate of the 
significance of the low rejection of the null result. 



9 SUMMARY AND CONCLUSION 

We have performed a cross-correlation analysis between the 
WMAP5 CMB data and various large-scale structure trac- 
ers including our new high redshift AAfi-LRG survey. The 
summarised conclusions of our findings are as follows: 

• We have found a null ISW result for z ^ 0.7 AAfi-LRG 
sample. The standard model is rejected at « 3% significance 
by this dataset. 

• We have confirmed the marginal correlations between 
WMAP5 CMB temperature fluctuations and SDSS LRGs 
at z=0.35 and 2SLAQ LRGs at z=0.55. 

• The null result in the AAJ7-LRG sample at large scales 
is unlikely to be caused by the negative contribution of the 
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SZ effect, given its angular extent and the expected ampli- 
tude of ISW signal. 

• We have made a range of tests on the AAfl cross- 
correlation measurement which confirms its robustness. 
These include moving the magnitude limits up to 0.5mag 
brighter, removing areas of sky with significant dust absorp- 
tion, using an estimate of the cross-correlation that takes 
out any possible systematic effects due to SDSS stripes and 
comparing the standard and uber-calibrations of the SDSS 
photometry. We have also checked the effects of stellar con- 
tamination in our samples. All these tests produce results 
consistent with our original measurements. 

• We have also reproduced the cross-correlation results 
of most previous authors using our techniques. In particu- 
lar we have reproduced the marginally positive correlations 
seen using SDSS magnitude limited samples of galaxies and 
NVSS radio sources. 

• However, rotation tests indicate that accidental align- 
ment or some unknown systematics can give rise to a cor- 
relation signal comparable to and in many cases even larger 
than the ISW signal itself. This suggests that the previous 
positive detections may still be subject to unknown system- 
atic effects. 

• Combining the new z « 0.7 LRG survey with the lower 
redshift LRG samples, the overall cross-correlation result 
is now as consistent with a null detection as it is with the 
standard ACDM model for both W-band and ILC data. 
For the ILC map, the significance of rejecting the standard 
model is « 2a whereas the result is only 0.5a away from the 
zero correlation hypothesis. 

• Given the results of the rotation test on the SDSS and 
2SLAQ LRG samples, the support these give to the stan- 
dard ISW model in the combined sample may have even 
less statistical weight than indicated above. 

• There is a possibility that the absence of the ISW cor- 
relation in the high redshift dataset is due to evolution of 
the dark energy equation of state. However, we regard it as 
unlikely that evolution could take place over the short red- 
shift interval between the 2SLAQ and AAfl datasets. It is 
more plausible that the differences between the redshift bins 
are purely statistical, particularly given the rotation test re- 
sults. We note that the individual positive detections that 
we have discussed are only marginally statistically signifi- 
cant and the combined ILC dataset is more consistent with 
zero than with the standard model prediction. 

• If the ISW effect was generally absent then the impact 
on cosmology would be large because this would be strong 
evidence against an accelerating Universe. This would there- 
fore argue against a significant role for a cosmological con- 
stant or dark energy in the Universe. Moreover, the absence 
of ISW would also argue against any modified gravity model 
which produced acceleration. The model which would be 
heavily favoured would be an Einstein-de Sitter model with 
O m = I. However, if such a model had a critical density of 
exotic, CDM particles then there might be a contradiction 
with the high baryon densities in rich galaxy clusters such 
as Coma. This rich cluster 'baryon catastrophe' has previ- 
ously argued against a high CDM density because starting 
from fij,/fl m ~ 0.03, it was difficult to understand in a hier- 
archical model how to produce a 5x big ger baryon fraction 
in rich galaxy clusters (|White et al.lll993T ). 

• It is therefore important to repeat the LRG measure- 



ments made here, now in the Southern sky. One opportunity 
to do this will arise from the new ESO imaging surveys in 
the South which are about to start, the VST ATLAS and the 
VISTA Hemisphere Survey. If the results we have found here 
are repeated then there could be significant consequences for 
cosmology. 
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